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\  primary  concern  of  Navy  management  is  the  ability  to  meet  the  Navy's 
manpower  requirements  in  terms  of  both  the  quantity  and  quality  of  personnel. 
Recruitment,  promotion,  retirement,  and  other  managerial  policies  are  all  directly 
related  to  and  restricted  by  the  Military  Personnel,  Navy  (MPN)  budget.  Because 
about  17  percent  of  tlie  MPN  budget  is  allocated  for  officer  basic  pay  (approximately 
$1  billion),  special  attention  must  be  given  to  the  accurate  forecasting  of  officer 
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Imputations.  The  Naval  Pay  Predictor,  Officer  (N/VPPO)  Model,  designed  to  arcom 
plish  this  objective,  relies  solely  upon  historical  population  data  and  user-supplied  end 
strengths.  Time  series  analysis  is  used  to  determine  a  general  set  of  forecasting 
models  that  adequately  explain  the  historical  data.  Other  statistical  procedures, 
including  those  employed  in  determining  the  cost  of  paying  the  officer  force  and  in 
estimating  newly  commissioned  officer  input  populations,  are  also  described.  Valida¬ 
tion  results  indicating  errors  of  less  than  .21  percent,  .64  percent,  and  1.06  percent  lor 
1-,  2-,  and  1-year  lead  times,  respectively,  lor  total  officer  basic  pay  are  presented. 
These  results  represent  an  improvement  of  more  titan  SO  percent  over  the  average  of 
the  last  9  years  of  forecasting  by  the  budget  planners. 
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FOREWORD 


1  he  effort  desi  ribed  in  tins  report  supports  the  development  of  compensation  and 
tost  models,  an  exploratory  development  objeetive  under  Task  Area  Zf: *>■>. 32 1 .0 10.  The 
ob|ei  tive  of  this  task  area  is  to  develop  techniques  to  forecast  the  cost  and  behavioi 
elfeits  of  the  Navy's  compensation  polities.  The  main  eflort  in  FY79  was  directed 
toward  solving  two  problems.  The  first  concerns  the  development  ol  quantitative  methods 
for  evaluating  alternative  retirement  systems.  The  objective  is  to  forecast  the  cost  and 
retention  behavior  effects  of  such  systems  on  the  enlisted  force.  The  second  problem  is 
forecasting  the  Navy's  military  manpower  budget  to  avoid  cost  overruns  and  personnel 
turbulence.  This  problem  has  been  addressed  by  developing  the  Naval  Personnel  Pay 
Piedictor,  Officer  (NAPPO)  model,  which  is  described  in  this  report,  and  the  companion 
Naval  Personnel  Pay  Predictor,  Enlisted  (NAPPE)  model  for  enlisted  personnel,  which  is 
.  urrently  being  used  by  the  Deputy  Chief  of  Naval  Operations  (MPA  T)  to  forecast  enlisted 
basic  pav  and  to  distribute  the  enlisted  force  over  3 1  length-of -service  categories.  The 
NAPPO  model  will  be  combined  with  the  NAl’Pl  model  and  other  forecasting  techniques 
to  yield  a  Military  Personnel,  Navy  (MPN)  account  forecasting  model,  to  be  used  by  OP- 
131,  the  Financial  Management  Aranch.  Liaison  will  be  maintained  with  OP-131  as  to 
when  the  all-MPN  forecasting  model  will  be  installed. 

This  report  is  the  second  in  a  series  documenting  work  in  the  area  of  budget 
forecasting.  The  first,  NPRDC  Technical  Report  78-4,  described  the  development  of  the 
NAPPE  model. 


DONALD  F.  PARKER 
Commanding  Officer 


SUMMARY 


Problem 

Navy  personnel  planning  and  policy  decisions  often  involve  problems  of  predicting 
future  events  or  conditions.  One  such  problem  is  that  of  forecasting  the  Navy’s  military 
manpower  budget  to  avoid  cost  overruns  and  personnel  policy  turbulence.  In  particular, 
cost  overruns  can  occur  due  to  inaccurate  predictions  of  either  the  personnel  force 
structure  (personnel  arrayed  by  length-of-service  (LOS)  and  pay  grade)  or  the  various  pay 
rates.  To  control  the  first  source  of  error,  forecasting  models  are  being  developed  to 
predict  required  obligations  based  on  a  desired  or  allocated  man-year  average  by  pay 
grade. 


ODiective 


The  objective  of  this  effort  was  to  provide  a  detailed  description  of  the  analysis  and 
procedures  used  to  formulate  a  model  for  forecasting  Navy  officer  force  structures  and 
basic  pay  obligations.  The  model,  which  will  be  used  to  monitor  the  service  age 
characteristics  of  the  force  and  its  costs  in  basic  pay,  is  known  as  the  Naval  Personnel 
Pay  Predictor,  Officer  (NAPPO).  — 


NAPPO  relies  upon  historical  USN,  USNR,  USNT,  and  All  Navy  (ALNAV)  quarterly 
force  structure  files  dating  back  to  1963.  Time  series  analysis  techniques  were  applied  to 
these  files  to  find  a  particular  set  of  time  series  models  that  would  be  appropriate  for 
forecasting  the  LOS  marginal  distribution  of  each  array.  Various  combinations  of  these 
forecasts  were  then  compared  to  obtain  a  "best"  forecast  for  the  ALNAV  LOS  distribu¬ 
tion.  Additional  statistical  procedures,  previously  developed  for  the  Naval  Personnel  Pay 
Predictor,  Enlisted  (NAPPE)  model,  were  used  for  (1)  deriving  the  interior  of  the  force 
structure  matrix  given  the  forecasted  LOS  and  inputted  pay  grade  marginal  distributions, 
(2)  forecasting  the  force  structure  for  personnel  with  less  than  1  year  of  service  or  more 
than  30  years,  (3)  costing  the  force  structure,  (<0  estimating  average  strength,  and  (5) 
validating  the  model. 

Findings  — - - 


The  statistical  techniques  employed  in  NAPPO  proved  to  be  highly  accurate  in 
producing  estimates  of  officer  basic  pay.  Validation  results  indicated  forecasting  errors 
of  less  than  .3  percent  for  FY76,  FY77,  and  FY78.  NAPPO's  predictions  for  mean  LOS  of 
the  force  also  indicated  a  high  degree  of  accuracy.  As  expected,  the  accuracy  of  the 
forecasts  generally  diminished  as  the  forecast  lead  time  increased.  A 


Conclusions 

1.  Sufficient  data  are  available  to  allow  time  series  analysis  techniques  to  be 
applied  toward  the  development  of  a  model  to  forecast  the  officer  force  structure  and, 
hence,  basic  pay  obligations. 

2.  A  particular  set  of  time  series  models  explain  a  great  majority  of  the  data.  This 
set  of  models  provides  a  solid  statistical  basis  for  making  force  structure  forecasts. 

3  Comparisons  show  that  the  NAPPO  model  is  much  more  accurate  than  the 
current  procedures  used  to  forecast  officer  basic  pay. 
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Recommendation 

The  NAPPO  model  should  be  used  by  Navy  management  to  provide  future  basic  pay 
and  mean  LOS  estimates. 
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INTRODUCTION 


Problem 

The  Navy's  ability  to  meet  its  manpower  requirements,  in  terms  of  both  quantity  and 
quality  of  personnel,  is  restricted  by  the  size  of  the  personnel  budget.  Recruitment, 
promotion,  retirement,  and  other  managerial  policies  are  all  directly  or  indirectly  related 
to  the  amount  of  dollars  available  in  the  Military  Personnel,  Navy  (MPN)  budget,  of  which 
17  percent,  or  approximately  1  billion  dollars,  is  allocated  for  basic  pay  of  officer 
personnel.  In  fiscal  year  1978,  $1,093,603,000  of  obligations  were  incurred  for  this  one 
item.  Thus,  even  a  1  percent  discrepancy  in  forecasting  the  officer  basic  pay  component 
may  result  in  an  over-  or  underexpenditure  of  $10  million.  Because  the  amount  of  basic 
pay  is  directly  proportional  to  the  officer  force  structure  (i.e.,  the  number  of  personnel 
tabulated  by  pay  grade  and  length  of  service),  the  force  structure's  configuration  must  be 
forecasted  accurately  to  predict  future  basic  pay  accurately. 

Background 

The  officer  force  structure  is  characterized  by  14  pay  grades  and  32  length-of- 
service  (LOS)  categories.  A  representative  force  structure  matrix  is  depicted  in  Table  1. 
The  pay  grades  are  composed  of  the  regular  officer  force  (0-1  through  0-10)  and  the  four 
warrant  officer  grades  (W-l  through  W-4).  The  LOS  category  refers  to  the  number  of 
years  in  service,  with  cell  0  including  all  personnel  with  less  than  1  year  of  service  and 
cell  31  including  all  personnel  with  31  or  more  years. 

Once  a  force  structure  matrix  is  projected,  a  simple  procedure  is  followed  to  obtain 
the  cost  associated  with  that  configuration.  First,  because  the  force  structure  fluctuates 
over  time,  the  projected  population  for  a  certain  period  of  time  is  obtained  by  finding  the 
average  number  of  personnel  within  a  pay  category  at  the  beginning  and  at  the  end  of  the 
period.  Next,  the  cost  is  calculated  by  multiplying  the  average  strength  within  each  pay 
category  by  a  prespecified  statutory  rate  and  summing  over  all  pay  categories.  The  125 
pay  categories  and  a  set  of  representative  rates  are  depicted  in  Table  2.  Not  included  in 
Table  2  are  the  pay  rates  for  0-1's,  0-2's,  and  0-3's  with  4  or  more  years  enlisted  service. 
Naval  Personnel  Pay  Predictor,  Officer  (NAPPO)  does  not  attempt  to  forecast  this 
particular  group  of  personnel  because  future  pay  grade  totals,  which  are  input  to  NAPPO, 
do  not  include  numbers  for  these  personnel.  The  primary  task,  then,  in  forecasting  the 
officer  basic  pay  entails  the  forecasting  of  the  force  structure. 

The  NAPPO  model  employs  several  sources  of  data.  First,  four  separate  quarterly 
force  structure  files  dating  back  to  1963  are  maintained.  The  four  files  contain  United  States 
Navy  (USN);  United  States  Navy,  Reserve  (USNR);  United  States  Navy,  Temporary 
(USNT);  and  total  force  (ALNAV)  data,  the  latter  being  the  sum  of  the  other  three  files 
for  corresponding  grades,  years  of  service,  and  point  in  time.  Thus,  snapshots  of  those 
forces  over  time  provide  the  basis  for  applying  time  series  analysis.  Also  used  by  NAPPO 
for  validation  and  forecasting  purposes  are  (1)  the  pay  tables  (e.g.,  Table  2)  that  have  been 
in  effect  since  1963,  (2)  managerial  inputs  in  the  form  of  pay  grade  totals  for  each  quarter 
that  is  to  be  predicted,  and  (3)  actual  (historical)  pay  grade  totals.  If  predictions  are  to  be 
reliable,  all  of  the  data  base  must  be  kept  as  up-to-date  as  possible. 

Objective 

The  objective  of  this  report  is  to  provide  a  detailed  description  of  the  analysis  and 
procedures  used  to  formulate  the  NAPPO  model. 
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APPROACH 


The  basic  approach  taken  by  NAPPO  is  an  aggregative  one.  Although  the  use  of  the 
historical  quarterly  inventories  allows  ready  observation  of  personnel  flows,  making  direct 
predictions  for  each  length  of  service  (LOS)  cell  would  require  additional  knowledge  of 
promotions,  such  as  tnne-in-grade  and  flow  points.  This  information,  in  turn,  would  not 
only  lead  to  the  expansion  of  the  data  base  but  also  to  a  much  more  complex  model. 
Instead,  NAPPO  relies  xdely  upon  historical  data  to  project  the  future.  It  forecasts  the 
LOS  marginal  distribution  by  considering  its  statistical  properties  and  then  distributes  the 
personnel  over  the  authorized  pay  grade  totals  to  arrive  at  the  final  force  structure 
matrix.  The  following  sections  describe  in  detail  the  methods  employed  to  forecast  future 
inventories. 

Forecasting  the  LOS  Distribution 

Examination  of  the  quarterly  force  structure  tables  suggests  two  alternative  ap¬ 
proaches  for  predicting  LOS  cells  1  through  31.  (Cell  0  is  a  special  case  and  is  discussed 
beginning  on  page  9.)  The  first  involves  the  use  of  the  actual  cell  populations,  denoted 
P.(t)  (i.e.,  the  population  with  longevity^ at  time  _t,  with  _t  being  measured  in  quarters).  As 

will  be  demonstrated  shortly,  a  forecast  for  LOS  cell  |at  time  _t,  F.(t),  is  calculated  using 

actual  and/or  previous  forecasts  for  the  quarters  through  time  t-1.  That  is, 


A 

F.(t)  =  P.(t-l). 


(1) 


The  second  approach  consists  of  using  a  set  of  transformed  variables  that  measure 
the  transition  of  the  force.  They  are  defined  as  net  loss  rates  (r)  and  are  given  by 


r.(t)  =  (P-(t)  -  Pj  +  1(um))/P.(t),  (2) 

for  j  =  0,1,.. .,30.  (The  loss  rate  for  cell  31  is  also  a  special  case  and  is  discussed  on  page 
11.)  Thus,  (l-r.(t)),  or  Cj(t),  is  the  continuance  or  transition  rate;  namely,  the  proportion 

of  personnel  in  LOS  cell  j_  at  time  X  that  "move"  to  cell  j  +  1  at  time  t  +  m.  Since  LOS  is 
measured  in  years,  m  is  set  equal  to  4  in  NAPPO,  thus  providing  yearly  transition  rates.  A 
forecast  for  LOS  cell  j  +  1  at  time  _t  is  then  computed  as 


Fj+1(t)  =  C.(t-4)P.(t-4)  (3) 

A 

where  C  .(t-4)  is  the  predicted  rate  of  flow  from  P^(t-4)  to  P- +  ^ (t).  Because  P.(t-4)  is 
known,  this  method  requires  only  the  forecasting  of  the  continuance  rates. 

Given  the  four  data  sources  available  in  the  USN,  USNR,  USNT,  and  ALNAV  files, 
there  are  four  straightforward  ways  of  forecasting  the  ALNAV  population  of  LOS  cell 
j  +  1  at  time  ;t  using  the  two  approaches  outlined  above.  They  can  be  written  as 


t-  Pj2(t-4)Cj2(t-4)  +  Pj3(t-4)Cj3(t-4), 
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Vw = y*-*  v-*1. 


Ff  ,(I|  >  R  ,  ,(t-l)  .  R  ,  ,<t-l)  .  R  ,  ,<t-l) 


or 


j+r  '  VI, i 


Fj\i(t)  -  WM) 


i+i.2' 


VI, 3' 


(4) 


where  P .  ^(t)  designates  the  actual  number  of  personnel  in  cell  j_  at  time  _t  for  data  set  k  = 

1,2, 3, 4  indicating  USN,  USNR,  USNT,  and  ALNAV,  respectively.  Note  that  the  first  two 
forecasts  use  predicted  continuance  rates,  whereas  the  latter  two  use  actual  population 
projections.  Also,  the  first  and  third  forecasts  combine  USN,  USNR,  and  USNT 
predictions  to  give  an  ALNAV  forecast. 


As  was  done  in  the  NAPPE  model,  time  series  analysis  techniques  (cf.,  Box  <5c  Jenkins, 
1970;  Brown,  1963)  were  applied  to  the  four  data  sets  to  forecast  the  continuance  rates. 
For  each  of  the  124  time  series  (31  different  loss  rates  for  each  of  the  four  data  files), 
each  consisting  of  58  data  points,  the  rates  were  plotted,  along  with  their  autocorrelation 
and  partial  autocorrelation  functions.  Also,  the  first  differences  of  the  loss  rates  and 
their  autocorrelation  and  partial  autocorrelation  functions  were  computed  and  plotted. 
Seasonality  factors  were  identified  by  calculating  the  autocorrelation  functions  and  their 
standard  errors.  Finally,  a  chi-square  test  was  made  on  the  first  24  autocorrelations  to 
determine  whether  they  could  be  distinguished  from  white  noise.  This  procedure  was  used 
in  an  attempt  to  identify  a  general  class  of  time  series  models  that  would  be  applicable  to 
most  or  all  of  those  of  the  individual  time  series  while,  at  the  same  time,  being 
conservative  in  the  consumption  of  computer  time. 


Table  3  presents  the  series'  seasonality  factors  and  suggests  one  or  more  appropriate 
time  series  models  based  on  the  autocorrelation  and  partial  autocorrelation  functions. 
The  terminology  employed  is  that  of  Box  and  Jenkins  ( 1 970,  Chapters  3  and  4),  where  a 
(p,d,q)  autoregressive  integrated  moving  average  (AR1MA)  model  is  differenced  d  times 
and  contains  p  autoregressive  and  q  moving  average  parameters.  Thus,  an  entry  in  the 
table  such  as  "(1,1,0)  x  4"  implies  that  a  (1,1,0)  model  be  fitted  to  the  series  created  by 

taking  Z  -  Zt  ^  for  t  >  5,  where  Zt  is  the  tth  observation  in  the  series.  This  step 

eliminates  a  seasonality  factor  of  period  4. 


The  undifferenced  data  in  a  large  majority  of  each  set  of  series  indicated  non¬ 
stationary  behavior.  Many  of  these  appear  to  be  well  suited  to  a  single  exponential 
smoothing  model  because,  in  general,  they  do  not  vary  about  a  fixed  mean,  yet  they 
exhibit  homogeneous  behavior  of  a  kind.  In  this  regard,  Box  and  Jenkins  (1970)  comment: 
".  .  .  although  the  general  level  about  which  fluctuations  are  occurring  may  be  different  at 
different  times,  the  broad  behavior  of  the  series,  when  differences  in  level  are  allowed 
for,  may  be  similar"  (p.  II).  The  exponential  smoothing  model  is  most  easily  recognized 
by  its  recursion  formula, 


Z{(a)  =  aZt  l  +  (l-a)Zt_ j(a),  for  0  <  a  <  1, 


(5) 
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*15  A  . 


which  shows  th.it  each  now  level  of  the  process  (the  forecast),  Z^lu),  is  arrived  at  hy 

interpolating  between  the  new  observation,  Z (  .,  and  the  previous  level,  Z}  j(ci).  With  a 

at  its  uppei  limit  ol  I,  the  level  ol  the  process  is  simply  the  latest  observation,  with  all 
previous  history  ignored.  With  or  small  (closer  to  0),  the  previous  level  is  given  more 
weight  than  the  latest  observation.  The  smaller  n  is,  the  longer  the  time  period  required 
by  the  model  to  recognize  a  <  hange  in  "levels"  ol  the  data. 

Seventy  seven  of  tire  time  set  irs  seem  to  l>e  modelled  well  hv  single  exponential 
smoothing.  First,  of  the'  senes  weie  reduced  to  white  noise  allot  first  differencing 
implying  a  (0,1,0)  model  (i.e.,  Z (  Z(  j  *  a(  where  a(  is  the  noise  component  and  Z^  is  the 

t*1’  observation  of  the  senes),  hut  this  is  simply  equation  ("*)  wither  equal  to  1.  As  shown 
m  Table  1,  these  series  wore  UxN  s «- 1 1  s  t  «>,  7-1  I,  IS  70,  77,  Z7,  2*1;  UNNR  cells  7,  6,  .  10. 

11,1  U 7 1 .  24,  25,  20,  10;  USNT  c  ells  10.  17  I  74.  25,  27,  2*h  '();  and  Al  NAV  c  ells  7.  1.  6. 
7,  I  1,  It,  20-22,  75-27.  2 It  was  indicated  that  another  70  senes  (I  isN  cells  14.  71; 
l.l.SNK  cells  *>,  22.  26.  77;  USNT  cells  I  6.  2<>;  and  ALNAV  cells  4,  10.  14.  16-18,  21)  could 
be  fitted  by  a  (0,1,1)  ARIMA  process  (i.e.,  Z  “-j  )  '  at  (1at  \'  vvtir,r‘  ^  1  Sm‘ *' 

this  is  the  inverted  form  ol  equation  (5)  (see  box  .V  Jenkins,  pp.  10"’  106),  howevei ,  those 
series  can  also  be  forecasted  using  the  single  exponential  smoothing  model.  Note  that 
double  and  triple  exponential  smoothing  are  not  applicable  here  Irecause  no  evidence 
exists  of  any  linear  oi  quadratic  trends  in  the  data. 

Another  '2  senes,  which  required  lust  differencing  to  achieve  stationai  ity,  contained 
evidence  of  a  seasonality  fac  tor.  A  significant  autocorrelation  at  lag  7  was  noticed  in  b 
of  these  series  (I  ISN  cell  28;  USNR  c  ell  \  USNT  cell  78:  and  Al  NAV  c  ells  0.  \  28).  Ifeie. 
a  significant  autocorrelation  at  lag  n  means  that  pairs  ol  observations  separated  by  n  time 
periods  are  significantly  correlated.  Twenty-six  series  possessed  an  indication  of  a  cycle 
of  length  4  (l  ISN  c  ells  11-17,  21,  25,  2<’.  10;  USNR  cells  1.  4.  8.  12.  21;  USNT  cells  1 1.  16- 
23;  and  ALNAV  c  ells  *>,  II.  12.  I*t.  10). 

Finally,  66  of  tlie  series  could  also  be  modelled  well  bv  an  autoregressive  model  of 
order  I,  the  (1,0.0)  model.  In  notational  form,  the  model  is  Z?  -  4>Zt  a^,  where  I  c|. 

I.  Table  )  indicates  the  <><•  sot  ies  to  be  USN  cells  4-6,  8  -21,  27,  7'h  USNR  c  ells  7,  5  IS, 

20,  22,  71,  25,  27,  28;  USNT  c  ells  0,  10-21;  and  At  NAV  cells  0.  7  4.  6,  *»,  II.  15.  Some  ol 
these  cells  are  such  that  two  or  more  time  series  models  were  suggested  by  the 
autocorrelation  functions.  This  was  the  case  with  manv  of  the  (1,0,0)  c  ells. 

because  the’  single  exponential  smoothing  model,  the  seasonal  model  of  period  4,  and 
flic  first  cv  dei  autoregressive  model  adequately  explain  76  of  the  USN  series,  7'*  of  the 
USNR  series,  2/  of  the  USNT  senes,  and  26  ol  the  Al  NAV  senes,  no  turther  models  foi 
the  loss  rates  were  considered.  The  remaining  18  series  not  covered  hv  these  models  are 
explained  by  an  assortment  of  various  models,  including  the  (0,0,1),  (2,0,0),  and  (1,1.0) 
ARIMA  processes.  Since  eac  h  of  tliese  models  requires  the  estimations  ot  at  least  one 
parameter,  increased  forecasting  accuracy  was  sacrificed  in  the'  inlet est  ol  using  less 
computer  time.  II  should  l>e  noted  that  the  a  parametei  in  the  exponential  smoothing 
model  is  estimated  bv  choosing  that  ct  ftom  the'  set  jo,  .05,  .10,  ...,  .**5,  I . ())  that 
minimizes  historical  error.  The  same  is  true  tot  the  parametei  in  the  single  ordei 
autoregressive  model  although,  technically  speaking,  could  xnrv  fiom  I  jo  I.  Negative 
values  fc’r  if  were  excluded  since  nc’  negative  autocorrelations  ot  lag  I  weie  ohsei  ved  in 
the  time  set  ies.  In  effect,  then,  21  exponential  smoothing  models  (tot  21  o's).  I 
autoregressive  models  (fot  21  ^'s),  and  I  seasonal  model  (requiring  nc’  parametei 
estimation)  are  fitted  to  all  124  time  sei  irs  and  the  model  that  best  explains  the  past  is 
chosen  to  forecast  future  transition  rates. 


,s 


, 


Similar  procedures  were  used  to  produce  a  set  of  models  that  would  fit  most  of  the 
actual  population  data.  Table  4  summarizes  the  analysis  for  the  124  time  series  covering 
LOS  cells  1-31.  (Cell  0  is  a  special  case  and  is  discussed  below).  In  contrast  to  the 
enlisted  LOS  cell  populations,  very  little  seasonality  is  evident  in  the  officer  populations. 
Ninety-two  of  the  124  time  series  are  best  modelled  by  the  single  exponential,  first  order 
autoregressive,  or  seasonal  model  of  period  4,  all  of  which  are  described  above.  The 
exponential  model  alone  appears  to  bo  adequate  for  79  of  those. 

USN,  USNR,  USNT,  and  ALNAV  forecasts  for  a  particular  quartet  are  provided  tor 
each  LOS  cell  using  continuance  rates  and  actual  population  data.  These  individual 
forecasts  are  produced  by  the  model  that  minimizes  historical  error.  These  eight 
forecasts  for  each  cell  (two  forecasts  for  a  particular  LOS  cell  in  each  ot  four  data  sets) 
are  used  to  produce  the  four  ALNAV  forecasts  seen  in  equation  (4),  where  the  USN, 
USNR,  and  USNT  forecasts  are  combined. 

Examination  of  the  forecasts  F1  through  I'1*  m  equation  (4)  showed  that  the  forecasts 
based  on  continuance  rates  (F1,  F‘)  yielded  less  historical  error  than  those  based  on  actual 
population  data  (FJ,  F").  Consequently,  in  some  cases,  the  latter  forecasts  could  be  used 
to  "damp"  the  continuance  rate  forecasts.  Such  instances  would  be  characterized  by 
conditions  in  which  tlie  continuance  rate  varies  inversely  with  the  source  population. 
Thus,  a  fifth  prediction  was  made  by  taking  a  weighted  sum  of  the  first  four: 


F5 

)+ 


,(t)  -  p(F*  ,(t)  r  F?  ,(t))/2 
I  F  1  U  l 


(1-B)  (F5  .(t)  ♦  F*  . ( t ))/ 2, 
I*  I  I* 1 


(6) 


0  <  6  <  l.  Here  B  is  some  weight  that  produces  a  minimum  sum  of  squared  error.  Finally, 
t lie  forecast  for  a  particular  cell  is  tlx'  one  chosen  from  { F 1 ,  F2,  F5,  I  1*.  !;S]  that  lias  been 
the  best  estimator  of  the  past,  with  "best"  implying  the  forecast  that  minimizes  historical 
error.  As  will  be  seen,  this  "final"  forecast  will  most  likely  be  modified. 

Special  Considerations  for  Cells  31  and  0 

As  noted  earlier,  LOS  cell  31  contains  all  personnel  with  more  than  3|  years  of 
service,  because  the  continuance  rate  as  previously  defined  only  explains  t he  movement 
of  personnel  from  cell  30  to  cell  31,  a  special  rate  was  defined  that  would  encompass  the 
movement  of  personnel  from  cell  30  as  well  as  the  continuance  of  those  personnel  already 
in  cell  31.  It  is  written 

r 3 1 ( 1 3  ,  ((P,0(t)  ♦  P31(t))  -  P3,(t*4))/(Pw(t)  ♦  P3,(t)).  (7) 


(Compare  with  (2).)  The  time  series  composed  of  these  rates  was  analyzed  in  the  same 
fashion  as  those  for  ^(t)  through  r^(t). 

Cell  0  is  also  a  special  case;  that  is,  no  continuance  rate  can  lx'  defined  for  it  because 
there  is  obviously  no  source  population  with  military  service  less  than  0  years.  Also,  due 
to  tlx'  fluctuating  input  of  newly  commissioned  officers  (determined  largely  by  personnel 
policies),  the  high  degree  of  instability  of  populations  in  coll  0  make  it  undesirable  to 
smooth  tlx-  historical  populations  in  that  cell.  Vet,  one  ready-made  forecast  is  available 
from  the  realization  that  the  sum  of  the  forecasts  for  cells  0  through  3|  must  equal  the 
total  force  size.  Thus,  we  have  the  "residual  forecast": 
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Seasonality  Factors  and  Appropriate  Models  for 
Actual  Population  Time  Series 


%  O  O  •"*  O  “  *  -  •  r-t  “  C  O  O  -  «  -r-4.~4r-4.-4r-4r-4.-4.-4 

H  !  •  •  •  •  O  •  *  iH  M  •  O  fH  H  r-4  •  r-4  •  •  •  O  O  O  •  -  •  •  «  m  m  m 

U&*  r~*  *-4  r~i  (D  “  r-4  r-4  -  *  O  •  •  •  »  H  -OOO  «  •  *r-4HrHHHHHM<N 

0.-0  -  -  O  O  -  r-4  O  O  O  •  O  •  »  .  r-4  r-4  r~4  .  . . 

O  O  C'  (N  O  h  — O  O  O  ■*—  v^v^o*-  fNrjrjw^v.  OOOOOOO  O  X 

fc*  SI  *»-'  x—  w  ^  ^  ^  •*  •  -  •»  •»  •*  •»  ^  •»  w  v--  •  •  •  •  •  •  v»-  w  n_-  »—  v/  v_-  v_. 

O.  O  O  ^  OOO  O  tN  ^4  ^  O 


UWliJblWUlUJUJUJUiUJ^UJUlUlUiUUJUlUlUl^UlUUbJblUUJUtN 

sc  •/.  z  ;*  z  y.  y.  z  z  v.  z  z  z  sc  y.  y.  y.  y.  y.  y.  y  y.  y.  y  y.  y.  y.  y.  y 

ooocooooooo  oooooooooooooooooo 

7.  z  z  y.  7.  7.  y  y  y.  y.  y.  y.  y.  y.  /.  y  y.  y.  y.  y.  z  y.  y  y.  y.  y  y  y  y 


0^0  -V  -  O  ^  -N  -  ^  — t  O  0  .~>  ,-V  ,~V  -X  .~S  - 

,  -  -4  -  O  <N  --400000^-4000000  -  -  -4  -ooooooo 

-*  -  #H  «  «•  — « . .  •  •  •  »QO  •  . 

-  O  -  O  O  -0—4— 4— 4— 4r-4r-4»— 4  OOOOOO  «  «  —4  -r_4r-4r-4tH.-4r-4.-4 

Jo  -  O  -  “  O  «  r-4  -•-••-• 

,‘-'0'-«'00^00000000<N<NICNCNI<NC'I^WO'—'— <»-«— 4—4—40.-4 


—4  u  I 

«fl  o  j  rs<  u:  UJ  uj  n:  r*  w  u>  u:  u;  u:  uj  ui  ui  ui  ui  u:  uJ  uj  u:  uj  u:  uj  u;  uj  ui  u:  uj  u>  w  uj  I 

c  ^  »*  z  «’  It.  y.  s 4  /.  z  sc  z  z  z  z  z  z  z  z  z  z  z  y.  z  z.  z  z  z  z  z.  z  , 

ou  oooo  ooooooooooooooooooooooooo 

2  *  ,  y  y.  7.  y.  y.  y.  y.  y.  y.  y  y  z  z  z  z  sc  sc  sc  y.  sc  sc  z  z  z  z  z  z  sc  z 

3  u- 


O  -«  -e-j^^oOOOOO^O 


41  H  H  H  -  H  HHHHHHHHHHhHH  »  H  »  -  »  *  -  »  H  H  H  H  H 

O.T?  »  -  «  O  . *  -  O  -000000 

O  o  OcmOv^OOCOOOOOOOOOOO'-^O'-^wwv^wn— 'OOOOO 

Li  JI  W  v—  •  «  W  V_wv_  ••  V^S^WSrV^S.  W  .  •  v_  •  »  ..  •  ■>  .  •  M.  .«  w  V_ 

O.  »-S  N  ,-S  N  4— S  S  4— \ 

CL  o  ooooooo 


OOOOOO 


UJ  uj  uJ  -r  uJ  ui  ui  u;  u)  b]  uj  u)  u.’  ui  to  to  u:  ui  u:  uj  b.  uJ  uj  uj  uj  io  to  ui  uj  bJ  to 

y.  y.  y.  y.  y.  y.  z.  z  sc  sc  y  y.  y.  y.  y.  y.  y.  z  z  z  z  Z  z  z  y.  y.  y.  y.  y.  y. 

OOO  OOO  300000000000000000000000 
z  z  z  z  z  z  sc  z  z  z  sc  z  z  z  z  z  z  z  z  z  Z  Z  Z  Z  Z  Z  7.  z  z  z 


Li  &J  -4<N4'"*-tr-4,_4— 4-« 

0.-0.  *  -  «  -  - 
OO  OXXXOOOO 

W.  s*  w 


O  O  •  H  »*«*•*  O  H  »  «  »  •  H  •  »  •  • 

•  *  O  •  H  H  H  rt  H  H  *  •  rH  H  r— 4  4—4  *  H  H  r- 4  H  tN 

r-4  f— 4  -O  *  »  -  “  »  »rHO  -«••••  o  •••• 

wwr-4^  000000v_v^0000v..  OOOO  x 


til  M  -f  ul  W  Kl  .1  vf  uj  Ul  UI  UJ  bl  111  UI  UJ  UJ  W  W  UJ  bl  W  W  U)  ui  U)  Ul  W  1.1  ' 

z  z  z  y  z  z  z  y.  z  z  z  z  z  z  z  z  z  z  z  z  z  z  z  z  z 

O  OOO  ■>  ooooooooooooooooooooo 

Z  Z  Z  z  z  z  sc  sc  z  z  z  s*.  sc  z  z  z  z  z  z  z  z  z  sc  z  sc 


10 


where  is  the  ALNAV  forecast  for  cell  0,  F  is  the  total  force,  and  P.  is  the  best  forecast 

selected  from  {F1,  F2,  F3,  F4,  I-'5}  for  L05  cell  j.  This  forecast  for  cell  0  may  be 
abnormally  small,  large,  or  even  negative,  however,  due  to  large  changes  in  the  total 


force  size  and  the  relative  stability  of  forecasts  through  P.^.  Therefore,  the  following 


analysis  was  performed  to  obtain  more  realistic  forecasts  for  all  the 


First,  proportions  were  computed  that  relate  the  size  of  cell  zero  to  the  total  force 
size  over  the  historical  data  base.  Symbolically,  these  ratios,  denoted  R.,  are  determined 
by  1 


Rj  =  P0(i)/F(i),  for  i  =  1,  ...,  62  (9) 

where  \_  represents  62  quarters  of  history  and  P^U)  and  F(i)  are  the  cell  0  and  total  force 

populations  for  quarter  L  Next,  the  R.'s  were  examined  for  trends  and  patterns  using  time 

series  analysis.  As  it  turned  out,  a  second-order  moving  average  process  applied  to  first 
differences  fits  the  data  extremely  well.  This  model  is  written  Z  =  Z  j  -<f>at  2,  where 

af  2>  t^ie  error  for  period  t-2,  is  the  actual  value  minus  the  forecasted  value  for  that 

period  and  <f>  lies  between  1  and  -1.  Thus,  the  forecast  for  a  cell  zero  to  total  force  ratio 
for  time  period  _t,  denoted  Rt,  is  determined  from 


Rt  -  Rf  |  -  .5666at  2 

or  (10) 

R,  =  Rt_j  -  .5666(Rt  2  -  Rt_2). 

A 

Here  d  .5666  and  R^  2  is  the  forecast  for  time  period  t-2.  Applying  the  forecasted  ratio 
to  the  user-supplied  total  officer  force  yields  another  forecast  for  the  cell  0  population.  In 
addition,  95  percent  confidence  limits  for  Rt  can  be  derived  by  substituting  the  95  percent 
confidence  limits  (.3^259,  .79061)  for  the  value  of  if>  in  equation  (10). 

The  final  formula  for  computing  the  cell  0  population  for  any  future  quarter  t  uses 
the  residual  forecast  for  that  quarter  as  well  as  equation  (10): 

.  +  (Rt  -  -5666a.  ,)*F(t)) 

P*(t)  =  max{(R t  -  .79061at_2)»F(t),  min  v  — - —  — - — - ’ 

(Rt  -  .34259at  2)*F(t)  )  }  (11) 


1 


where  the  first  term  in  the  brackets  is  the  lower  bound  on  the  size  of  cell  0;  the  second, 
tlx*  average  of  the  residual  forecast  and  the  moving  average  forecast;  and  the  third,  the 
upper  bound  on  the  size  of  cell  0.  The  lower  and  upper  bounds  in  equation  (11)  assume 
a^_2  >  0.  If  a{  2  '•  *he  coefficients  .34239  and  .79061  would  be  interchanged. 


It  is  quite  likely,  however,  that  the  adjusted  cell  0  lorecast,  P*(t),  is  different  from 

^(t),  t lx*  residual  forecast.  Thus,  the  discrepancy,  l^(t)  -  P*(t),  must  be  distributed  over 

each  of  tlx'  remaining  cells.  The  procedure  used  ensures  that  the  cells  that  have  been 
predicted  with  a  high  degree  of  accuracy  are  altered  only  slightly.  It  distributes  the 
discrepancy  according  to  the  amount  of  the  prediction  error  in  a  cell  relative  to  the  total 
prediction  error  over  all  other  cells: 


1! 


PjUt) 


iv(c)  -  (i; 


t)  -  P^(t))*(cj(t)Pj(t)/ 


i 


(12) 


j  -  1,2,. ..,3 1,  where  the  error  is  calculated  as 


t-l 

Ej(t)  =  ^|Pj4(i)-P)(i>|/Pj4(.)  (13) 

for  cell  i .  Recall  that  P)4(i)  denotes  the  actual  ALNAV  population  for  cell  ^  in  quarter 

Also  note  that  four  quarters  are  required  to  initiate  the  forecasting  process  (for 
computation  of  continuance  rates);  hence,  the  lower  range  in  the  summation. 

Having  executed  tlx*  above  procedure,  NAPPO  now  has  the  LOS  marginal  distribution 

(Pg.  Pj.  P  >»  •••.  p  ^ | '  oi  the  force  structure  matrix.  The  vector  is  complete  and  sums  to 

the  total  authorized  strength.  The  pay  grade  vector  is  either  supplied  by  the  user  or 
interpolated  by  NAPPO  (discussed  on  page  14).  Thus,  the  problem  now  is  to  calculate  the 
interior  entries  in  the  force  structure  matrix.  Once  this  has  been  accomplished, 
multiquarter  forecasts  are  obtained  by  concatenating  the  predicted  quarter  with  the 
historical  data  base  and  employing  the  optimal  time  series  models  that  were  chosen  for 
the  first  prediction. 

Constructing  the  Force  Structure  Matrix 

Given  the  two  marginal  distributions  of  the  force  structure  matrix,  NAPPO  follows 
•Hosteller's  (1968)  procedure  for  deriving  the  interior.  This  iterative  procedure  involves 
the  renormalization  of  the  rows  and  columns  of  a  standard  matrix  to  obtain  the  given 
marginals,  while  maintaining  the  relative  "associations"  among  the  entries  of  the  standard 
matrix.  In  NAPPO,  the  standard  matrix  is  computed  as  the  average  of  the  previous  12 
historical  ALNAV  inventories.  Although  the  Hosteller  procedure  is  documented  in  the 
NAPPE  report  (Chiptnan,  1977),  it  is  presented  here  again  for  convenience. 

To  begin,  suppose  that  the  LOS  and  pay  grade  marginals  are  denoted  by  L;  and  G  , 

respectively,  where  i  =  0,1,.. .,31  and  j  i  1,2,..., 14.  With  these  and  the  standard  matrix 
being  given,  Hosteller's  procedure  involves  the  following  steps: 


12 


1.  Let  n  =  I. 


2.  Let  S"*1  -  (Lj/ES")  S?  (or  all  u 
,  I  ‘I 

\  Let  s"(*2  ,  (Gj/ES^*  ’>S^  1  for  all  j. 


4.  If  T.'. S  k  L  for  all  i  go  to  (2),  otherwise  stop. 

)  '' 


This  process  allows  estimates  of  a  force  structure  matrix  to  be  made  for  specific  points  in 
time.  To  predict  basic  pay  costs  for  a  tiscal  year,  some  estimate  of  tlie  configuration  of 
that  matrix  throughout  the  year  is  necessary. 


Estima 1 1 ng  Average  Strength 


Given  an  actual  inventory  at  some  quarterly  interval  and  a  forecast  inventory 
representing  a  period  four  quarters  later,  an  estimate  of  average  strength  can  be 
computed;  that  is, 


SM  .  (SM  .  V/2. 


where  5^  and  represent  the  actual  and  forecasted  strength  matrix  entries  respectively. 
Clearly,  the  average  force  5  derived  in  this  manner  assumes  that  changes  in  a  force 
configuration  over  four  quarters  occur  in  a  constant  relation  to  time.  Thus,  if  the 

A 

difference  between  T.Z  S  and  rZ  S  were  100,  an  average  matrix  would  reflect  a  monthly 

lj  ij  1 

change  of  8.33  or  a  quarterly  change  of  23,  etc. 

Straightiine  methods  are  inadequate,  however,  because  changes  do  not  occur 
uniformly  throughout  the  year.  Unfortunately,  the  only  practical  alternatives  irTterms  of 
available  data  are  average  strength  computations  based  on  quarterly  increments.  Thus, 

A 

given  S()  for  t  r  December  1973,  a  forecast  of  S.(  is  made  for  t  =  December  1976. 

Similarly,  four-quarter  forecasts  are  made  Tune  to  June,  September  to  September,  and 
March  to  March.  So,  by  forming  force  structure  matrices  for  each  quarterly  prediction 
and  by  combining  those  four  forecasts  with  the  most  recent  actual  matrix,  a  relatively 
good  estimate  of  a  yearly  average  strength  can  be  computed  by  a  five-point  average: 


A  A 


Up  to  this  point,  estimates  of  future  pay  grade  distributions  were  assumed  to  be  an 
exogeneous  input  accessible  to  the  users  of  the  model.  The  estimation  of  the  14  cells  of 
the  pay  grade  marginal  distribution,  however,  is  no  trivial  task. 


While  the  end-year  pay  grade  structure  is  roughly  established  by  the  budget 
(authorized  end  strength),  the  quarterly  pay  grade  distributions  are  largely  a  function  of 
management  actions  taken  throughout  the  fiscal  year.  If  pay  grade  totals  for  September, 
December,  March,  and  June  are  input  by  the  user,  the  model  will  form  complete  force 
structure  matrices  based  on  those  given  pay  grade  marginals  and  the  forecasted  longevity 
marginals.  This  input  will  be  based  on  official  plans  that  reflect  management  intentions. 
Specifically,  the  pay  grade  populations  for  the  0-3  through  0-10  and  W-2  through  W-4 
mainly  reflect  the  planned  phasing  of  promotions  during  the  subject  year,  while  0-1,  0-2, 
and  W-l  primarily  reflect  the  planned  phasing  of  input. 

If,  instead,  the  pay  grade  totals  for  the  three  quarters  interior  to  the  fiscal  year  (the 
beginning  total  is  known  and  the  ending  total  is  always  input  as  a  matter  of 
strength/budget  policy)  are  to  be  forecasted  by  NAPPO,  the  following  procedures,  akin  to 

Aw  /\ 

those  in  NAPPE,  are  used:  Let  F  .  and  F  .  .  be  the  actual  or  user  input  pay  grade  totals 

for  the  end  (or  beginning)  of  2  consecutive  fiscal  years.  Estimates  of  the  quarterly  pay 
grade  totals  are  then  given  by 


and 


/V  A  A 


't.l,i  =  (Ft,i 


A  A  A 


(16) 


These  simple  linear  interpolations  could  be  altered  by  increasing  or  decreasing  them  from 
the  straightline  results  by  reference  to  an  average  of  historical  seasonal  patterns  of  pay 
grade  population  movement.  Since  analysis  revealed  wide  variability  in  the  quarterly 
deviations  from  straightline  averages  from  year  to  year,  however,  it  was  decided  to 
allow  the  user  to  implement  any  seasonal  pattern  deemed  appropriate  by  directly 
inputting  the  quarterly  future  pay  grade  totals.  Otherwise,  NAPPO  will  use  those 
obtained  in  equation  (16). 

Costing  the  Force  Structure 

The  procedures  outlined  previously  produce  a  forecast  of  the  force  structure  matrix 
for  the  end  of  each  quarter  (31  December,  31  March,  30  June,  and  30  September)  during 
the  forecast  year.  In  addition,  the  beginning  matrix  (30  September)  is  available  in  the 
form  of  known  (i.e.,  historical)  data. 

The  following  costing  procedures  are  observed.  If  the  pay  scale  is  constant  for  the 
year,  the  five  matrices  or  inventories  can  be  averaged  and  a  pay  scale  such  as  that  in 
Table  2  applied  to  the  average  force  to  produce  a  forecast  of  annual  basic  pay  costs.  If 
the  pay  scale  changes  during  a  fiscal  year,  some  provision  must  be  made  for  costing  part 
of  the  average  strength  with  one  pay  scale  and  for  costing  its  complement  with  another. 
Since  no  data  are  available  to  forecast  population  changes  within  a  quarter,  it  is  assumed 
that  intraquarter  changes  oc<  ur  uniformly  throughout  the  quarter.  Also,  if  the  pay  scale 
is  in  effect  for  a  full  quarter,  the  average  strength  for  that  quarter  is  defined  by  the 
average  of  the  begin  and  end  quarter  matrices,  and  the  cost  is  derived  by  applying  the 
constant  pay  rates  to  the  quarterly  average.  Finally,  if  the  pay  scale  changes  during  the 
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quarter,  then  (I)  the  force  structure  at  the  point  of  change  is  approximated  by  linear 
interpolation  between  the  begin  and  end  quarter  matrices,  (2)  the  populations  for  the  start 
and  end  of  each  pay  scale  are  averaged  to  obtain  the  average  force  to  which  the  relevant 
pay  scale  is  applied,  and  (3)  a  corresponding  daily  pay  scale  is  applied  to  the  average 
force.  The  results  then  are  multiplied  by  the  number  of  days  for  which  the  pay  scale  is 
applicable.  The  results  for  all  quarters—whether  single  or  dual  pay  scales—are  summed 
over  the  fiscal  year  to  produce  an  estimate  of  annual  basic  pay. 

This  method  is  followed  for  any  fiscal  year  to  be  predicted.  It  requires  only  the  end 
strength  pay  grade  totals  for  each  year  to  be  forecasted.  The  user  may  supply  quarterly 
pay  grade  totals  as  desired. 
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RESULTS 


Historical  Validation  of  the  NAPPO  Model 

One  problem  in  obtaining  a  complete  validation  of  NAPPO  lies  in  the  fact  that  the 
authorized  end  strengths  are  seldom  achieved  exactly.  Consequently,  although  some 
discrepancy  may  exist  between  the  actual  and  the  forecasted  force  structures,  it  would 
not  significantly  affect  the  validation  results. 

To  demonstrate  NAPPO's  accuracy,  the  model  was  used  to  forecast  historical  force 
structures  and  their  respective  costs  for  the  fiscal  years  1974-1978.  NAPPO  was  allowed 
to  use  historical  data  up  to  1,  2,  and  3  years  preceding  the  year  to  be  predicted.  All  pay 
grade  totals  used  were  actual  quarterly  totals  as  opposed  to  authorized  end  strength 
totals,  which  may  or  may  not  have  been  achieved.  It  should  be  noted  that,  although  the 
forecast  for  any  given  cell  in  the  force  structure  matrix  could  be  erroneous,  the  total  cost 
may  be  accurate  due  to  offsetting  errors  in  other  cells. 

Therefore,  an  additional  validation  criterion  was  established  to  test  for  systematic 
bias  in  the  forecasting  p.  ocess.  This  involved  comparing  the  forecasted  and  actual 
longevity  distributions,  with  the  criterion  being  average  length  of  service.  Table  3 
summarizes  these  results,  given  the  three  different  lead  times.  As  shown,  there  appears 
to  be  no  evidence  of  consistent  under-  or  overestimation  of  mean  LOS,  although,  in 
general,  the  forecasts  are  less  accurate  as  lead  time  increases. 

The  historical  validation  for  basic  pay  involved  the  comparison  of  the  costs  of  the 
projected  and  the  actual  force  structures.  The  absolute  values  of  the  percentage 
differences  between  the  two  also  are  included  in  Table  5,  which  demonstrates  NAPPO's 
accuracy  for  forecasts  made  1,  2,  and  3  years  in  advance  of  the  projected  period.  Overall, 
the  average  absolute  1-year  lead  time  fiscal  year  percent  error  in  cost  prediction  in  Table 
5  is  .21  percent,  with  a  standard  deviation  of  .68  percent.  Similar  figures  for  the  2-year 
and  3-year  lead  time  forecasts  are  .64  (1.55)  and  1.06  (1.96)  respectively.  In  forecasting 
mean  LOS,  the  average  fiscal  year  absolute  error  for  the  1-,  2-,  and  3-year  lead  times  are 
.09,  .25,  and  .45  years  respectively.  Once  again,  for  both  cost  and  mean  LOS  forecasts, 
accuracy  generally  diminishes  with  an  increase  in  lead  time. 

NAPPO  vs.  Current  Procedures 

One  means  of  measuring  the  usefulness  of  NAPPO  is  to  compare  the  NAPPO 
forecasts  with  those  in  the  annual  budget  submissions.  Each  year  the  budget  submission 
gives  I-  and  2-year  forecasts  for  all  of  the  Military  Personnel,  Navy  budget  activities, 
including  officer  pays  and  allowances. 

Figure  1  depicts  the  percentage  differences  of  the  NAPPO  and  budget  submission  2- 
year  forecasts  from  the  actual  obligations  incurred.  In  all  cases  NAPPO  was  more 
accurate,  with  an  average  absolute  error  of  .64  percent,  as  compared  to  4.00  percent  for 
the  budget  submissions.  The  1-year  forecast  errors  for  the  same  years  were  .18  percent 
for  NAPPO  and  1.40  percent  for  the  submissions.  When  the  data  for  FY7 1  and  FY72  (the 
poor  forecasting  years  for  the  submissions)  are  excluded,  NAPPO  is  still,  on  the  average, 
3-1/2  times  more  accurate  than  the  submission  in  2-year  forecasts  of  officer  basic  pay 
(.53%  vs.  1.96%).  In  FY78  dollars,  the  average  errors  represent  differences  of  $5.8  million 
and  $21.6  million  respectively. 
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Table  5  (Continued) 


CT' 

CN 

o 

sp 

cn 

HJ 

O 

PH 

vX? 

\X> 

CN 

CN 

O 

<r 

r-«. 

in 

in 

oo 

NX? 

00 

PH 

cn 

m 

PH 

<r 

in 

PH 

PH 

PH 

PH 

<y 

PH 

PH 

CN 

PH 

m 

in 

<r 

in 

in 

vX> 

00 

vX> 

vO 

sD 

O' 

VX> 

3  u 

H  W  U)  sU 
O  ^ 
•o  K  u  cu 


vC  00  (N  (M  rs 
O  w-4  <T  CN 


ph  ph  ao  co  O 

sj  ^  iTi  ID  K) 


CT>  in  O  rH  H 

CN  n  O  (N  H 


iJH  rH  O  CN  rH  rH 
(N  H  O'  (D  m  CN 


rH  O  0C  <r 
C  00  fD  (D  H 


+  +  +  I  + 


TJ 

•  0)  cn 

cn  *j  o 

^  u  J 
><  «H 
TO  C 
ph  a)  nj 
u  d) 
a.  T 


vj  ^  vf  O' 
H  O'  C  O'  J' 


vj  H  o  oo  vj 
^  >J  >J  vj  <T 


00  C*H  rH  CN  PJ 

in  m  uh  in 


h  rv  rH  vi'  rs 

CN  rH  CM  O  O  I-H 


r-  ao  V£>  O 

r-H  rH  C*  O 


rH  O  rH  O  O 


H  W  Cfi  0) 

o  ^ 


r-  cn  >£>  m 

(N  lD  N  \C  ID 


n  ao  (N  vO 
rH  O  CN  rH  o 


CN  00  ^  H  in 
m  i/n  rH 


O  VT*  O  ON  rH 

(N  fN  n  cn  cn  O 


pn  rH  p-.  ph 
H  H  (N  (N  O 


D  K  U  d, 
CO 


+  +  I  •  I 


I  I  I  rH 

I 


rs  >N  CN  v£*  (N 
m  cn  ph  m  rn 


rH  r^.  in  m  ph 
VT  CD  (*H  CD  nJ 


ph  in  rH  m 
CN  H  vj  O  CN 


(D  o  ID  m  (N 
CO  rH  CN  O  O  f-H 


n  vT  in  sc 

v£>  vO  v£>  \C  v£> 


W  cfl 
o  U 
KUO, 
T> 

<0 

i 

-4  T3 
0)  OH 
•  *J  O 
U  U  *-J 

>M  «rH 

"O  C 

f-H  0)  CT3 

U  0) 


r —  PH  ID  rH 
rH  O  »-H  O  C 


rH  CN  PH  sT  VD 
O  f-H  f-H  PH  fH 


•o  cc  O'  (n  n 
t-H  CN  rH  O  rH 


O'  00  PH  o  O'  CD 
O  O  rH  CN  rH  O 


ID  N  J>  vj  O' 
O  O  CN  r-  PH 


+  +  +  I 


•  +  +  I 


+  +  +  +  + 


r^vjO-sjcor^ 
m  <r  m  vj  -vj 


r-i  O  vO 
iD  sj  vO  CD  vj 


cn  O'  O  O 

Sf  -O  VO  ID  ST 


CD  vj  \C  vj  CN  -sf 
sj  \T  iD  vO  1^  iT 


<r  ph  m  o  o 
i —  r-^  m  ph  o 


ph  cn  >j  <r  m  m  m  d  in  iO  \C  iO  vT  it  i —  n  co  oo  oo  cc 


I 


11.09 


The  Implementation  of  NAPPO 

NAPPO  is  currently  operating  in  a  demand  mode  on  a  UNIVAC  I  I  10  computer  located 
at  the  Naval  Ocean  Systems  Center,  San  Hiego.  Currently,  plans  call  for  NAPPO  to  be 
included,  along  with  NAPPE,  in  a  comprehensive  budget  forecasting  model.  This  model 
will  provide  forecasts  for  all  budget  activities  of  the  Military  Personnel,  Navy  account, 
including  items  such  as  enlisted  and  officer  basic  pay,  enlisted  and  officer  allowances, 
subsistence,  and  permanent  change  of  station  pay.  A  future  report  will  document  the  all- 
MPN  model,  particularly  the  techniques  used  to  forecast  the  budget  items  other  than 
basic  pay. 

As  a  stand-alone  model,  NAPPO  is  useful  not  only  for  forecasting  basic  pay  but  also 
for  monitoring  the  age  of  the  force  through  the  mean  LOS  forecasts.  Thus,  Navy 
management  can  be  kept  aware  of  upward  or  downward  trends  in  the  experience  of  the 
force.  Also,  for  general  use,  NAPPO  subroutines  allow  the  user  to  forecast  LOS  transition 
rates,  to  update  the  historical  data  base,  and  to  validate  the  model  with  respect  to  data 
ranging  back  to  19fO.  Still  other  routines  permit  the  examination  of  historical  transition 
rates  or  the  force  structure  matrices  themselves. 
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CONCLUSIONS 


1.  Sufficient  data  are  available  to  allow  time  series  analysis  techniques  to  be 
applied  toward  the  development  of  a  model  to  forecast  the  officer  force  structure  and, 
thus,  basic  pay  obligations. 

2.  A  particular  set  of  time  series  models  exists  that  adequately  explains  a  great 
inajority  of  the  data.  This  set  of  models  provides  a  solid  statistical  basis  for  making  forc  e 
structure  forecasts. 

3.  Comparisons  show  that  the  NAPPO  model  is  more  accurate  than  the  current 
procedures  used  to  forecast  officer  basic  pay. 

RECOMMENDATION 

It  is  recommended  that  the  model  be  used  by  Navy  management  to  provide  future  basic 
pay  and  mean  LOS  estimates. 
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